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Abstract. An integrable semi-discretization of the Camassa-Holm equation is presented. 
The keys of its construction are bilinear forms and determinant structure of solutions of the 
CH equation. Determinant formulas of A^-soliton solutions of the continuous and semi-discrete 
Camassa-Holm equations are presented. Based on determinant formulas, we can generate 
multi-soliton, multi-cuspon and multi-soliton-cuspon solutions. Numerical computations 
using the integrable semi-discrete Camassa-Holm equation are performed. It is shown that the 
■ integrable semi-discrete Camassa-Holm equation gives very accurate numerical results even 

in the cases of cuspon-cuspon and soliton-cuspon interactions. The numerical computation for 
an initial value condition, which is not an exact solution, is also presented. 
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1. Introduction 

The Camassa-Holm (CH) equation 

wt + 2k wx — wtxx + 3wwx = Iwxwxx + wwxxx ■ ( 1 • 1 ) 

has attracted considerable interest since it has been derived as a model equation for shallow- 
water waves [jll. Here, w = w{X, T), K is a positive parameter and the subscripts T and X 
appended to T) denote partial differentiation. Originally, this equation has been found in 
a mathematical search of recursion operators connected with the integrable partial differential 
equations [|2]. The CH equation has been shown to be completely integrable. In the case of 
K = 0, the CH equation admits peakon solutions which are represented by piecewise analytic 
functions [3]. Schiff obtained single- and two-soliton solutions in a parametric form by using 
the Backlund transformation [4]. An approach based on the inverse scattering transform 
method (1ST) provides an explicit form of the inverse mapping in terms of Wronskian 
ttHmiTllHl. The A^-soliton solution was also constructed by using the Hirota bilinear method 
[El [ini [m [13 (See also [[Bl). The key point of computations of the A^-soliton solution by 
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the Hirota bilinear method in those papers is the relationship between the CH equation and 
the AKNS shallow water wave equation [T4]. When K ^ 0, cusped solitary wave solutions, 
as well as analytic soli ton solutions, were found in [|l?l[T6l[T7l . In [[T8l[T9l . the interaction of 
cusped soliton (cuspon) was studied in detail. 

It is extremely difficult to perform numerical computations of the CH equation due to the 
singularities of cuspon and peakon solutions. So far, several numerical computations of the 
CH equation were presented [|20l |2ll l22l |23l IM 1251 |26l. However, none of these numerical 
methods gives satisfactory results for soliton-cuspon and cuspon-cuspon interactions. 

Integrable discretizations of soliton equations have received considerable attention 
recently [|27ll28l . Ablowitz and Ladik proposed how to construct integrable discrete analogues 
of soliton equations based on Lax pairs Il29ll30l . Hirota proposed another method to construct 
integrable discrete analogues of soliton equations based on bilinear equations [3T1 |32l l33l . 
Applications of integrable discretizations of soliton equations were considered in various 
fields [I34l,l35l.f36l. 

The purpose in the present paper is to present an integrable semi-discretization of the 
CH equation through a bilinear approach and to perform numerical computations by using 
integrable scheme. We show that the integrable scheme gives very accurate numerical results 
even for the soliton- solition and cuspon-cuspon interactions. 

The outline of this paper is as follows. In §2, we give bilinear forms and a determinant 
formula of the //-soliton solution of equation (|l.ll) . As far as we are concerned, bilinear 
forms directly related to the CH equation have not been known yet. These bilinear forms 
can help us to understand mathematical structure of the CH equation more thoroughly. Based 
on these bilinear equations, we obtain the determinant formula for the CH equation by using 
determinant technique [|37l l38l [39l l40l . In §3, we give an integrable semi-discrete Camassa- 
Holm equation and its determinant formula of the A/^-soliton solution. In §4, we present results 
of numerical computations by using the proposed integrable semi-discrete CH equation. 

2. Bilinear equations and determinant formulas of the Camassa-Holm equation 

In this section, we give bilinear equations and the A^-soliton solution of the CH equation. All of 
the existing works in constructing A^-soliton solutions of the CH equation using Hirota bilinear 
method take advantage of the relationship between the CH equation and the AKNS shallow 
water wave equation [I9l[10l[lilll2l[l4l. So far, the bilinear forms which are directly related to 
the CH equation remain unknown. In the present paper, we give bilinear equations obtained 
from the CH equation directly and derive the determinant solution by using the determinant 
technique. Our formulation in this section is crucial for the derivation of integrable semi- 
discretization of the CH equation. 
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Lemma 2.1. Bilinear equations 

^^D,D,-l^f-f^-gh, 
^Dt{Dy-2cD^)f-f=-D^g-h, 
^D,D,ff = D,gh, 



- 2cD,) -2^f-f= {DtD, -2)g-h, 



(2.1) 



have a determinant solution 



where 



In = 



= Pix + pjy + —^t + - — 2 5 + ^iO , 
Pi-c [Pi-cY 



(n) 

in) (w+1) 
¥2 ¥2 ' 



^ ¥r^^ - ¥r^-^^ 



2 1 

y\i = qiX + q^y^ 1 + 



qt - c (qi - c) 
Proof. Consider the following Casorati determinant solution, 

M (n+l) (n+N-l) 

¥r ¥1 ■•■ ¥1 ' 

¥2 ¥2 ' ¥2 



2^ + ^10- 



in) (n+l) 

¥V Tn ^ 



(n+iV-l) 



(2.2) 



where \|/- '^'s are arbitrary functions of four continuous independent variables, x, y, t and 5, 
which satisfy linear dispersion relations. 



(«+2) , ^ (n+l) , 2 (n) 



(n) 

Thus we can choose \|/- as follows 



(2.3) 

(2.4) 
(2.5) 
(2.6) 
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= piX + p]y H —t + - — 5 + ^,0 , 

Pi-c [pi-cY 

2 1 1 

V[i = qix + -J H 1 + + ri,o . 

qt-c [qi-cY 

For simplicity, we introduce a convenient notation, 



|ni,n2,---,"A'| 



(2.7) 



In this notation, the solution for the above bilinear forms, x„, is rewritten as 

%n = \n,n + \,---,n + N-\\. (2.8) 

We show that the above x„ actually satisfies the bilinear equations (12.11) by using the Laplace 
expansion technique Il39ll40ll . 

The differential formulas for x are given by 

{d^-Nc)Xn = \n,n + \,---,n + N-2,n+N\, (2.9) 

{dy-2cd^ + Nc^)%n = \n,n+\,--- ,n+N -2,n + N +l\ 
- \n,n+l,---,n + N-3,n+N- l,n + N\, 



(2.10) 
(2.11) 
(2.12) 
(2.13) 

(2.14) 



dtXn = |n — l,n + 1, . . ■ ,n+N — 1 1 , 

dsXn = \n — 2,n + 1, - ■ ■ ,n + N — 1\ — \n — l,n,n + 2, - ■ ■ ,n + N — 1\ , 
{dt{dj,~Nc)-l)Xn = \n-\,n + \,---,n+N-2,n+N\, 

dt{dy-2cdx+Nc^)%n = \n-l,n+l,---,n + N-2,n + N+l\ 

- \n-l,n + l,---,n+N-3,n + N-l,n+N\ , 

ds{dx-Nc)Xn = \n-2,n+l,---,n + N-2,n + N\ 

- \n-\,n,n + 2,--- ,n + N-2,n + N\, (2.15) 

{ds{dy-2cdjc+Nc^)-2)Xn = \n-2,n+l,---,n+N-2,n + N+l\ 

-\n-l,n,n + 2,--- ,n + N-2,n + N+l\- \n-2,n + \,- ■ ■ ,n + N -3,n+N - l,n + A^| 
+ \n-l,n,n + 2,---,n + N-3,n + N-l,n+N\ , (2.16) 

which are proved by using the linear dispersion relations (I2.3I) - (I2.6I) . (See Appendix) 



The first equation of eqs. (12.11) 

Let us introduce an identity for 2N x 2N determinant, 

n-l i n + l ... n+N-2 i n+N \n\ i n+N-l 

n-l i i n+N i n i n + \ ■■■ n+N-2 i n+N-l 



0. 
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Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear identity 
for determinants, 

|/i-l,n + l,---,n+A^-2,/7 + A^| x + 1, • ■ ■ ,n + A^- 2,n +A^- 1 1 
-\n,n+\, - ■ ■ ,n + N -2,n-\-N\ x |n - l,n + 1, • ■ ■ ,n + A^- 2,n +Af - 1 1 
+ |/7 + l,---,n+A^-2,n + A^-l,/7+A^| x |/7 - l,n,n + 1, ■ ■ +A^- 2| =0, (2.17) 

which is rewritten by using (12.81) . (|2.9I) . (|2.1 II) and (12.131) . into the differential bilinear 
equation, 

{dt{d^-Nc) - 1)X„ X X„ - {d^-Nc)Xn X 3?X„ + X„+iX„_i =0, 

i.e., 

(3f3xX„-x„)x,i -3fX,i3xX„ + x„+ix„_i =0. 

Setting n = 0, / = Xo ,g = Xi ,/z = X-i , the above bilinear equation leads to the first equation 
of (1211). 

The second equation of eqs. (12.11) 

Let us introduce two identities for IN x IN determinants. 



n-l I n + 1 ■■■ n+N-3 n + N-1 \ n+N \ n \ 

n-1 I I n+N \ n \ n + l ■■■ n + N-l 

n-l \ n + l ■■■ n+N-2 \ n+N+l \ \ n + N-l 



0, 



0. 



n-l \ \ n+N+l \ n n + l ■■■ n + N-l \ n+N-l 

Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear 
identities for determinants, 

\n-l,n + l,---,n + N -3,n+N - l,n + N\)< \n,n + l,- ■ ■ ,n + N -2,n+N - l\ 

- \n,n + l,---,n+N -3,n+N - l,n-fA^| x \n-l,n+l,- ■ ■ ,n+N -2,n + N -l\ 
+ \n-l,n,n + l,--- ,n + N -li^n+N -l \ x \n+l,n + 2,- ■ ■ ,n+N -l,n + N\ =0, 

(2.18) 

\n-l,n + l,---,n + N -2,n+N +l \ x \n,n + l,-- ■ ,n + N -2,n+N -l\ 

- \n,n + l,--- ,n+N -2,n+N +l \ x \n-l,n+l, - ■ ■ ,n+N -2,n + N -l\ 

+ \n+l,--- ,n + N -l,n+N +l \ x \n-l,n,n+l,- ■ ■ ,n+N -?,,n + N -2\ =0. 

(2.19) 

Taking the difference of these two bilinear identities, it is rewritten by using (|2.8I) - (|2.1 II) and 
(|2.14l) into the differential bilinear equation, 

dt{dy - 2cd^ + Nc^)Xn X X„ - [dy - 2cdx + Nc^)Xn x 3,x„ 

+ {d^-Nc)%n+l ^^n-l-^n+\{^x-Nc)Xn-l =0, 
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i.e., 

{dtdyXn)Xn - dyXn^t%i " 2c( (3,3xX„)x„ - 3;,X„3fX„) + (3xX„+l)x„_i - X„+i [d/Zn-l] 



which is nothing but the second equation of (1271 



The third equation of eqs. (12.11) 

Let us introduce two identities for IN x IN determinants. 



n — 2 n + 1 



n + N-2 i n+N i n 



I n+N \ n \ n + l 

n-l n n + 2 ■■■ n+N-2 \ n+N \ n + l 



i n+N-l 

' =0, 

■■■ n + N-2 \ n + N- 1 

i 



I n+A^ I n + l \ n n + 2 ■■■ n + N-2 \ n+N-l 

Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear 
identities for determinants, 

\n-2,n + l,n + 2,---,n+N-2,n + N\ x \n,n + 1,- ■ ■ ,n + N -2,n+N - 1\ 

- \n-2,n+l,n + 2,---,n+N-2,n + N-l \ x \n,n + l,n + 2,- ■ ■ ,n + N -2,n+N\ 

+ \n-2,n,n + l,n + 2,---,n+N-2\ x \n + l,n + 2,- ■ ■ ,n+N -2,n + N - l,n+N\ =0, 

(2.20) 

\n-l,n,n + 2,---,n + N-2,n+N\ x \n,n + 1,- ■ ■ ,n+N -2,n + N - 1\ 

- \n-l,n,n + 2,---,n + N-2,n+N-l \ x \n,n + l,n + 2,- ■ ■ ,n+N -2,n + N\ 

+ \n-l,n,n + l,n + 2,---,n+N-2\ x \n,n + 2,- ■ ■ ,n + N -2,n+N - l,n + A^| =0. 



0. 



(2.21) 



Taking the difference of these two bilinear identities, it is rewritten by using (|2.12l) and (|2.15l) 
into the differential bilinear equation, 

3^(3x-A^c)x„ xx„-3,-x„ X (3x-A^c)x„ + (3fX„_i)x„+i -x„_i(3rX„+i) =0, 



I.e., 



(3Axn)x„ -3i.x„3xX„ - (3;X„+i)x„_i -l-x„+i(3fX„_i) =0. 



Setting n = 0, / = Xo ,^ = Xi ,/z = x_i , the above bilinear equation leads to the third equation 
indU])- 



The fourth equation of eqs. (12.11) 

Let us introduce four identities for 2N x 2N determinants, 
n-2 i n + l ■■■ n+N-2 i n+N+1 \n\ 

n-2 i i n+N+1 i n i n + l ■■■ n+N-2 n + N-l 
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n—2 I n+l 



n + N-3 n+N-l I n + N I n 







n-2 
n—\\n n+2 



i n+N \ n \ n + l ■■■ n+N~ 1 

n + N-2 i n+N+l i n+l i 



0, 



n-l I 

n — 1 i n n + 2 



I n+N+l I n+l I n n + 2 ■■■ n + N-l 

■■ n+N-3 n+N-l i n+N i n + l i 



n — 1 i i n+N i n + l i n n + 2 ■■■ n 

Applying the Laplace expansion to the left-hand side, we obtain the algebraic 
identities for determinants, 

|n-2,n + I,- - ■,n + A^-2,n+A^ + 1| x |n,n + 1, ■ • • ,n + A^- 2,n +A^ - 1 1 
n,n + I,- - ■,n+A^-2,n+A^+ 1| x |n - 2,n + 1, ■ ■ ■ ,n +A^- 2,n + Af - 1 1 
n-2,n,n + l,---,n + A^-2|x |n + 1, ■ ■ ■ ,n + A^- 2,n +A^- l,n + A^+l| =0, 



+ A^-1 
bilinear 



0. 



+ 



+ 



+ 



+ 



(2.22) 



n-2,n + l,---,n + A^-3,n+A^- l,n + A^| x |n,n + 1, ■ ■ ■ ,n + A^- 1 1 
n,n + l,---,n+A^-3,n+A^-l,n + A^| x |n - 2,n + 1, • • ■ ,n +A^- 1 1 
n-2,n,n + l,---,n + A^-3,n+A^-l| x |n + 1, ■ ■ ■ ,n +A^- l,n +A^| =0, 

n-l,n,n + 2,---,n + A^-2,n+A^+l| x |n,n+l,---,n + A^-l| 
n,n + l,---,n+A^-2,n+A^+l| x |n-l,n,n + 2,---,n + A^-l| 
n-l,n,n + l,---,n + A^-2| x |n,n + 2,---,n+A^-l,n + A^+l| =0, 

n-l, n, n + 2, n + N-3, n+N-l, n+N\ x |n,n+l,---,n+A^-l| 
n,n + l,---,n+A^-3,n+A^-l,n+A^| x |n-l,n,n + 2,---,n + A^-l| 
n-l,n,n + l,---,n + A^-3,n+A^-l| x |n,n + 2, ■ ■ ■ ,n + A^- l,n +A^| = 0. 



(2.23) 



(2.24) 



(2.25) 



Taking an appropriate linear combination of these four bilinear identities, it is rewritten by 
using (12.161 ) into the differential bilinear equation, 

{ds{dy - 2cdx + Nc^) - 2)x„ X x,, - {dy - 2cd^+Nc^)Xn X d^Xn + 'dt^n-l X (3;, - A^c)x„+i 

- {dt{d^-Nc) - l)Xn-l X X„+l -Xn-l{M^x-Nc) - 1)X„+1 + (a^-A^c)X„_l X dfXn+l =0, 

i.e., 

- ((^f(^xX„+i)x„_i +3.^x„+i3fX„_i +3fX„+i3;,x„_i -x„+i3f3;cX„_i +2x„+ix„_i = 0, 

which leads to the fourth equation in (12.11 ). □ 
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Theorem 2.2. Bilinear equations 

-(^-D,D,-\^f-f^gh, 
2cff^{D, + 2c)g-h, 
[ -2ff={DtD, + 2cDt-2)g-h, 
have a determinant solution 



where 



f = ^0: 



h = X- 



(n) (n+l) (n+A^-1) 



1 



\i^PiXAr^—-t 
1 



r]i=-pix- 



-t+' 



1 



Pj + c' ' + 

Proof. In the previous Lemma, apply the 2-reduction condition 



This condition gives a constraint 



into bilinear equations. Thus we have 



1 



= -Pi ■ 



Dy = 0, 



DtD,-l]f-f=-gh, 



cDtDJf = D,gh, 
U)sDJ-f = Drg-h, 



^ {-cD,D,-2)f-f^{DtD,-2)g-h. 

After simple manipulations, we have 

1 



-D,D,-ljf-f^gh, 
2cff={D, + 2c)g-h, 
[ -2ff={DtD^ + 2cDt-2)g-h. 
Let us consider a determinant solution. When we apply the 2-reduction condition 
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on the determinant solution in the previous Lemma, we will have 

(") ("+1) (n+N-i) 

¥i ¥i ••• ¥i 

in) (n+\) (n+N-i) 
¥2 V2 ¥2 



(n) (n+l) 



where 



t,i = PiX H 1 + + ^,0 , 

Pi-c [pi-cy 
1 1 

V[i = -piX —t + - — —^s + ri,o . 

Pi + c {Pi + cy 

Thus the theorem was proved. 
Theorem 2.3. The CH equation 

(dj + wdx) {wxx -w)+ 2wx {wxx - w - K^) = , 

i.e., 

Wt + 2k wx — wtxx + 3wwx = '^wxwxx + wwxxx, 
where = l/c, is decomposed into bilinear equations 

' -(^-D,D,-l^f-f = gh, 

2cff={D, + 2c)g-h, 
[ -2ff={DrD, + 2cDr-2)g-h, 
through the hodograph transformation 

X = 2cx + logy, 
h 

T = t, 

and the dependent variable transformation 

Proof. Consider the dependent variable transformation 



u = — . 
f 



h 



From bilinear equations (|2.29l) . we obtain 

-((log/)xf-l) = wv, 
2c = {Dx + 2c)u • V, 

-2={DtD, + {2\ogfyt + 2cDt-2)u-v. 



(2.29) 



(2.30) 
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Eliminating (log/)^? from the third equation using the first equation in (12.301 ). we obtain 

2c = {Dx-\-2c)u-v , 

(2.31) 

—2 = [DtDx + 2cDt — 2uv) u-v. 
Thus we have 

2c = ( (log - ) + 2c ) uv , 

\\ V / X y ^2 32) 

-2 = (^{\oguv)j,t + (^log ^ j (log " j + 2c (log ^ j - 2uv^ uv . 
Eliminating (log + 2c from the second equation using the first equation in (12.321) . we have 



2c=((log'^ +2c^uv, 

-2 = ( (logMv);cf + — flog-) -2uv\uv. 
\ uv \ vJ t J 



Introducing new dependent variables 



^ S ^ gh 
V h 



we have 



2c 

— = (log(^)x + 2c, 

p (2.33) 
-2 = p(logp);,f + 2c(log^), -2p2. 



Let 



w= (log4))f = (log|)^ . 
From the first equation in (12.331) . we obtain 



Thus we have 



-(lOgp)f = ^W;,, 

;^(logp).,f+w+- = -p2. 
2c c c 



(2.34) 



Replacing (logp)r in the second equation by —^Wx using the first equation, we obtain 

-(logp)f = ^Wx, 

P P \ 1 1 2 

2c\Y^^) +''+c = cP • 

ZC V ZC ' X c c 



Consider the hodograph transformation 



X = 2cx + log(^, 
T = t. 
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Then we have 



11 



f ^ 2. + (l„g«, ^ 2c ^ (log l)-2cL^'j. (log«, ^ (log f )^ ^ w, 



5x = , 
P 

dt = dT + wdx ■ 

Using these results, we obtain 

-(3r + w3x)logp = wx, 

1 1 2 

-WXX + W+ - = -p • 
c c 

Eliminating p from the first equation using the second equation in eqs. (l2.35l) . we obtain 

(3r +w8x)log ( -w;cx +W+ - ] =-lwx- 



(2.35) 



Thus we finally obtain the CH equation 



(3r + wdx) {wxx - w) + Iwx ( wxx ~ ^ ~ ~ ) = ^ • 



1 



□ 



The theorem was approved. 

Furthermore we have the following theorem. 
Corollary 2.4. The CH equation has a determinant form. ofN-soliton solutions. 
Proof. From Theorem 2.2 and 2.3, the proof is obvious. □ 



3. A semi-discrete Camassa-Holm equation 

Lemma 3.1. Bilinear equations 

-Dt-l\f{k+l,l)- f{k, l) = -g{k+l,l)h{k,l) 
\ 

-A - 1 fik, ^ + 1) ■ m I) = -8{k, I + l)Kk, I) 




-D, -Dr\f{k+l,l)- f{k, l)=Drg{k+l,l)-h{k,l), 

-Ds - Dt \ f{k, / + 1 ) ■ fik, l)=Dtg{kJ + l)-h{k,l), 
have a determinant solution 

f{k, I) = Xo(fc, /) , g{K I) = xi (K I) , h{k, I) = x_i (fc, Z) , 



(3.1) 
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12 



where 



y^f\k,l) = ai,i{pi-cY{l-api)-\l-bpi)-'e^'+ai,2{qi-cf{l-aqiy\\-bqi)-'e^-^ ^ 



1 1 £ 



Pi-c {Pi-cf 
Proof. Consider the following Casorati determinant solution, 

^^^\k,i) y^^t'\k,i) ... ^vt^'-'Hkj) 



Xn{k,l) 



(3.4) 
(3.5) 



where \|/- ' 's are arbitrary functions of two continuous independent variables, t and s, and two 
discrete independent variables, k and /, which satisfy the linear dispersion relations, 

A,x^W=x^(''+^)+cvf\ (3.2) 

A,xk!"^=V^i)+cv|"\ (3.3) 

(n) (d-l) 

where and A; are defined as Ak^VikJ) = ^^^'^^"J^^"^'^^ and A/X|/(fc,/) = "Vik'O-jikJ-i) ^ 

(n) 

respectively. Thus we can choose \|/- ^ as follows: 

^l"\k,l) = {pi-cr{l-api)-\l-bpiy'e^' + {qi-cni-aqi)''{l-bqir'e'^', 

el Is 

Pi-c {pi-cy 
1 1 

Tl« = 1 + "2 5 + ri,-o . 



- c (qi - c) 
We use the following notation for simplicity 



X|/};^H/:i,/i) ¥!;'^(^2,/2) 



.(«2), 
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In this notation, x„(^, /) is rewritten as 

or suppressing the index k and /, 

Xn{kJ) = + 1, ■ --^n+N- 1| . 

The first equation of eqs. (13.11) 

By using the Casoratian technique developed in [|37l and [|38l . it is possible to derive the 
following differential and difference formulas for the x function, 

x„_i(fc,/) = \n-l,n,---,n+N-2\, (3.6) 
Xn{k+\J) = _^^^^_^ \nk+un+ \, - ■ ■ ,n + N -2,n+N - \k+i\, 

^ -\n,---,n+N-2,n + N-lk+i\, (3.7) 



(1-ac)^- 

(1 — ac) 



Xn+i{k+lJ) = _ ^ ^ \n + l,---,n+N- l,n + Nk+i\ 



^ 'n+l,---,n + N-l,n+N-lk+i\, (3.8) 



a(l-ac)^-2 

dtXn= \n-l,n+l,---,n+N-2,n + N-l\, (3.9) 

^—^Xn{k+lJ) = —. — ^—T^\n,---,n+N-2,n + N-lk+i\, (3.10) 
a a{l—acy^ ^ 

I — ac 



^ -dtXn{k+l,l)~rn{k+l,l) 
I — ac . 

= ~T^ rrr-j n- U.+i,/7 + l, ■■■,«+ A^- 2, n + A^- U+1 

ayV—acy^ ^ 

- Ife+i, ■■■,« + A^-2a:+i,«+ A^- 1/t+il 

^ |n- l,n+ l,---,n+A^-2,n + A^- l;t+l| 



a(l-ac)^-2' 
1 

(1 — ac) 



1 

l,---,n + A^-2,n+A^- 



(1-ac)^ 

= ^1 ^T^|n-l,n+l,---,n+A^-2,n + A^-U+i|, (3.11) 

a(l — ac)^^~^ 

Let us introduce an identity for 2A^ x 2A^ determinant, 

n+U ■■• n + N-2k \ n+N-lk \ n+N-lk+i \ n-lk \ nk \ 

I n+A^-l^ I n+A^-l^+i ; n-1^ ; n^ \ n + l^ ■■■ n+N-2j, 

Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear identity 
for determinants, 

|«-l,n+l,---,n + A^-2,n+A^-lyt+i| x \n,n + ■ ■ ,n + N -2,n + N - l\ 



0. 
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-|n,n+l,---,n + A^-2,n+A^-l;t+l| x |n-l,n+l,---,n+A^-2,n+A^-l| 

+ \n + \,---,n+N -2,n + N -\,n+N -\t-^i\ x \n-\,n,n + \, - ■ ■ ,n + N -2\ =0, (3.12) 

which is rewritten into the differential bilinear equation, 

—dtXn{k+\,l) -X„(fc+1,/) ) Xn{k,l)-- —Xn{k+\,l)dtXnik,l)+Xn+\{k+\,l)Xn-\{k,l) = . 



Setting n = 0, / = To , ^ = Xi , = x_i , the above bilinear equation leads to the first equation 
indTTl). 

The second equation of eqs. (13.11) 

The proof is similar to the proof of the first equation. 

The third equation of eqs. (13.11) 

We use the following differential and difference formulas for the x function 

x„_i(fc,/) = |n- l,n, ■■■,« + A^-2|, (3.13) 
x„(fc+l,Z) = \n,,+i,n+h+i,-- ■ ,n^N -2k+un^N -Ik+il (3.14) 
X„+i(fc+l,/) = |n + lA:+i,---,«+A^-U+l,'^ + A^fe+l|, 

= -^\n + \k+x,---,n+N - \k+i,n+N- 1|, (3.15) 
dtXn{kJ) = \n-\,n + \,---,n+N-2,n + N-l\, (3.16) 

- — —Xn{k+l,l) = -\n,n + \k+i, - ■ ■ -2k+un+N - lk+i\ 

a a 

= -\nk+un+l,n + 2k+u---,n + N-2k+un+N-\k+l\, (3.17) 
a 

dtTnik+lJ) = \n-lk+un+lk+u---,n + N~2k+un+N-lk+i\, (3.18) 

dtXn+l{k+ 1,Z) = \nk+un + 2k+u---,n + N-\k+un+Nk+i\ (3.19) 

= — \nk+un + 2k+u---,n + N-\k+un+N-l\, (3.20) 
a 

dtXn-\ik,l) = \n-2,n,---,n + N-3,n + N-2\, (3.21) 

dsXn{k,l) = \n-2,n + l,---,n + N-2,n+N-l\ 

+ \n,n-\,n + 2,---,n+N-2,n + N-\\, (3.22) 
1 — ac 1 — ac 

d,Xnik+ 1,/) = \n-2k+i,n + lk+ir ■ ■ ,n+N -2k+i,n + N - lk+i\ 

a a 

1 — ac 

H \nk+un- \k+un + 2k+u - ■ ■ ^n + N -2k+un+N - U+i , 

a 

= \n-lk+un + lk+ir--,n+N-2k+un + N-lk+i\ 
+ -\n-2,n + lk+u---,n+N-2k+un+N- 

--\n-l,nk+i,n + 2k+i,---,n+N-2k+i,n + N-lk+i\, (3.23) 
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Let us introduce following identities for IN x IN determinant, 

n-\k \ rik+i n + 2k+i ■■■ n+N-\k+i I «i I w+U i i n+N-\k 

\ rik \ n+lk \ n + 2k ■■■ n+N -lu \ n+N-lt 



n-lk\ 
and 

n-2k\ n + l^+i 



■■■ n+A^-l^+i i nk \ i n+A^-1^ 

i nk \ n + h ■■■ n+N-2k i n+N-lk 



0. 



n-2k 



Applying the Laplace expansion to the left-hand side, we obtain the algebraic bilinear 
identities for determinants, we obtain the algebraic bilinear identities for determinants, 

\n- l,nyt+i,n-|-2^+i,---,n + A^- U+i] x \n,n + l,n + 2, - ■ ■ ,n+N -2,n+N - l\ 
-\n,nk+i,n + 2k+i, - ■■ ,n+N - x \n- l,n + I,- ■ ■ ,n+N -2,n + N - l\ 
+ \n + l,nk+i,n + 2k+i,---jn + N- U+i] x \n~ l,n,n + 2, - ■ ■ ,n+N -2,n+N - l\ 
-|/iyt+i,n + 2yt+i, ■■■,n+N- Ik+i^n + N - 1| x \n - \ ,n,n + l,n + 2,- ■ ■ ,n + N -2\ 

= 0, 



(3.24) 



and 

\n — 2,n + Ij^j^i,-- ■ ,n + N —lk+i \ x \n,n+l,n + 2,---,n + N — 2,n+N—l\ 
— \n,n + - ■ ■ ,n+N — \k+i \ x \n — 2,n + \,---,n+N — 2,n + N—\\ 
+ \n + lk+U---,n + N~ lk+un+N-l\ x \n~2,n,n + I,- ■ ■ ,n + N -2\ 

= 0, 



(3.25) 



By dividing the difference of above two equations by a, we arrive at a bilinear equation 

1 — ac 1 — ac 

3,x„(fc+ \,l)Xn{k,l) x„(fc+ l,Z)a,x„(fc,Z) 

a a 

-dtXn{k+l,l)Xn{kJ)+Xn{k+lJ)dtXn{kJ) 

- dtXn+l (k+l, l)Xn-\{k, /) + X„+l (fc + 1 , l)dtXn-l{k, I) = . 

Setting n = 0, / = Xo , ^ = Xi , /? = X 1 , the above bilinear equation leads to the third equation 
in (1311). 



The fourth equation of eqs. O.ID 

The proof is similar to the proof of the third equation. 

Theorem 3.2. Bilinear equations 



I — ac 

Dt 

a 




fk+l ■ fk = - 


-gk+lhk, 


\ +ac 
a 




fk+l ■ fk = - 


-gkh+i , 


\ — ac 

D, 

a 


-Dt 


1 fk+l ■ fk = 


Dfgk+i-h 


l+ac 

Ds 


+ Dt 


fk+l ■ fk = 


Dtgk ■ hk+l 



(3.26) 
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have a determinant solution 



f{k, I) = Xo{k, 0) , g{k, I) = xi {k, 0) , hik, I) = x-i{k, 0) 

>+yv-l) , 



T„(A:,/) = 



(n+N-l) 
Y2 



(kj) 
{k,l) 



yV^\k,l) ^v\,^^'\k,l) 



where 

\|/|"^(^,/) = ai^i{pi-c)''{l-api)-\l+api)-^e^'+ai^2{-Pi-c)"{l+api)-\l -ap 



1 1 



Pi-c {pi-cY 



1 1 

t+- — ^2*+'n'0- 



Pi + c ipi^-c) 
Proof. Applying the 2-reduction condition 

qi = -Pi, b^-a, 

the T function satisfies 



1 



N 



-Xn{k,l). 



Y[{i-api){l +api) 
i=l 



Let 



Then 



fk = f{k,0), gk = g(k,0), hk = h{k,0), 
f \—ac \ 

\ — : — A - 1 1 A+1 • fk = -gk+ih , 
-A - 1 fk+i ■ fk = -gkh+i , 



a 

l+ac . 



a 

1 — ac 
a 

l+ac 
a 



-Ds - A j fk+i ■ fk = D,gk+i ■ hk , 
A + A ) fk+i ■ fk = Dtgk -hk+i. 



Let us consider a determinant solution. When we apply the 2-reduction condition 

<li = -Pi, b=-a, 
on the determinant solution in the previous Lemma, we will have 



x„{k,l) 



(n+l). 



>+i) 



(n+iV-1) 



(n+N-l) 
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where 

x\f\"\kJ)^ai,i{pi-cr{l-api)-''{l+api)-'e^'+ai,2{-Pi-cy\l+api)-'{l-api)-^e'^^, 



1 1 

-t + - 



Pi-c {Pi-cf^ 

1 1 
Pi + c {pi + cY^ 

Thus the theorem was proved. 

We propose a semi-discrete analogue of the CH equation 

K2K45?-4a2' 



dth= (l-^ ]h^wk, 



where a difference operator A and an average operator M are defined as 



k+l 



MFk 



Fk+i+Fk 



Theorem 3.3. The semi-discrete CH equation 



is decomposed into bilinear equations 
^ ^ \—ac 



5t K4-a2 



a 

\ +ac 



a 



l—ac 
a 

1 +ac 
a 



A - 1 1 fk+i ■ fk = -gk+ih , 
A - 1 j fk+i ■ fk = -gkh+i , 
Ds-Dtj fk+i ■ fk = Dtgk+i ■ hk , 
A + A j fk+\ ■ fk = Dtgk ■ hk+i , 



through the transformation 



4a 



(K2 + a)|^ + (K2-a)g: 



+1 
fk 



and 



where — Xjc. 



Wk = log 



gk 



hkJt 



(3.27) 



(3.28) 
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Proof. Let us start from bilinear equations (I3.28I ). By simple manipulation of eqs. (13.281 ). we 
obtain 

-2 - 1^ A+i ■ A = gk+ih+gkh+i , 
2acfk+\fk = (1 +ac)gk+\hk-{ \ -ac)gkhk+i , 
-2afk+ifk = ((1 + ac)Dt - a)gk+i ■ hk- {{I - ac)Dt + a)gk -hk+i. 



Let 



Then we have 



gk 
fk 



Vk 



h 
fk 



2 - log 



k+i 
fk 



1 = Uk+lVk + UkVk+l 



lac = (1 + ac)uk+iVk — (1 — ac)ukVk+\ 



-2a 



(1 +ac) 
(1 — ac) 



log 
log 



Hk+l 

Vk 
Uk 

Vk+\ 



+ log 



- log 



fk+\ 

fk 
fk+l 

fk 



— a 



+ a 



Uk+lVk 



^kVk+l ■ 



Substituting the first equation into the third equation in eqs. (|3.29l) . we obtain 

2ac = {1 + ac)uk+iVk - {1 - ac)ukVk+i , 



-la = (1 +ac)uk+iVk ( log 



Vk 



(1 -ac)ukVk+i ( log 



Uk 
Vk+l 



-:^iuk+iVk + UkVk+ \ + ac)uk+iVk+ {I - ac) UkVk+ 1 ) 
+ a^c{uk+iVk-UkVk+i). 

Simplifying the above equations, we have 

2ac = {l+ac)uk+iVk-{l -ac)ukVk+\ , 

-2a = ]-{{ \ + ac)uk+\Vk + [l - ac)ukVk+\) ( log "^"^^^^^^ ) 

^ V ^kVk Jt 



+ ac log 



Uk+\Uk_ 
Vk+lVk 



Let 



^k 



ifk_gk 
Vk h ' 

From the first equation in (13.301) . we obtain 



-2auk+\Vk+\UkVk: 
gkhk 



2ac ^k+i 
= ( 1 + j — ( 1 — ac) , 

UkVk+l ^k 

Multiplying these two equations, we obtain 



Pk = UkVk 



lac 



Uk+lVk 



fl 



(1 +ac) — (1 — ac^ 



'yfc+l 



(2ac)2 
p/t+ip/t 



(1 +ac 



(1-ac) ) ( (l+ac)-(l-ac' 



^k^ 



(3.29) 



(3.30) 



An integrable semi-discretization of the Camassa-Holm equation 
From the second equation of (13.301 ). we obtain 



19 



-2a = ]-{{\ + ac)uk+\Vk + {I - ac)ukVk+\) (log^^ ) + ac (log -2ap^,+ipA:. 
2 \ 9k J t 



Thus we have 

( 1 + ac) 1 Vyfc + ( 1 - ac) UkVk+ 1 
Aac 

Let us define a lattice parameter 

5^ 



log 



Pfc+i 



2 c c 



4ac 



(1 + ac)M^,+ iVA: + (l -ac)M;tV/t+i ' 



Then 



= 2 



Thus we have 



^ (1 +ac)Mfc+iV;^ - {\-ac)ukVk+i _ + ac)gk+\hk- { \ - ac)gkhk+i 
'(1 +ac)uk+iVk + {\ -ac)ukVk+i {\ + ac)gk+\hk + {I - ac)gkhk+i 

l+ac^k+\ ^ 
[l + ac)(Sfk+\ ~ - ac)(Sfk ^ ^ 1 - 
(l+ac)4)fc+i + (l-ac)4)^ ~ 1+ac 4)^+1 ^ ^ " 

1 - ac ^k 



l+ac 

2 



where a lattice parameter 6yt is a function depending on {k,t). The lattice parameter 
dX 2c 

corresponds to — = — = 2c+ (log(|))^ in the continuous case. At time t,X =Xq+ V 5yt is 

p 

a jc-coordinate of the k-th lattice point. Thus we have the following system: 
f (2ac)2 



Pk+ipk 



( 1 + ac) — ( 1 — '^c) ) ( ( 1 + ac) — ( 1 — ac 



Pfc , . 
1 +ac (^k+i 



+ o(log4'A:+l4»/t)f + - = -pyt+lP/t 

2 c c 



(3.31) 



J, l-ac (i}k 
1+ac 

1 - ac d^k 



-1 



+ 1 



Let 



Wk = (log(|);t)f = log 



^8k 



From the first equation in (13.311) . we obtain 

1 +flC / ^k+\ 

-(logpfc+lp^), = 



1 — ac \ <^k 



l+ac ^k+i 



<^k+ij, 



l-ac ^k 



-1 



1+ac (Sfk 
l-ac d^k+i 
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( l+ac(|);t+i \ 



I —ac (^k 



■ + 



l+ac(|)yt+i ^ 1+ac 



\l — ac (^k 
l+ac ^k+i 



loj 



l-ac (^k+ij 



1— ac (pyt 
l+ac(^k+]__^ 

l—ac (|)yt 



log 



^k+l 
<^k 



From the third equation in (13.311) . we obtain 

1+ac <^k+i 



2__^ ^ ^ l-ac 



l+ac ^k+\ 



I 



2 



I +ac ^k+i 



l—ac 

Multiplying these two equations, we obtain 

1-1= * 

8 



l—ac 



-I 



I +ac^k+i ^\ L 1 -ac (Sfk 



l-ac (Sfk 
Thus we have a system 



1 +ac ^k+i 



A{1 - ac){l + ac) 
{lacf 
9k+i9k 



-(logp;t+ipyt)/ = 2 



Wk+\ -Wk 



Let 



Then we have 



^ log-— + ^ + - = -pfc+ipfc 

Ok\ 9k J t 2 c c 



gkhk 



r^. = logp^ = log— 



-5f(''/t+l+'"yt) =2 



.2^2 



^Of('"/t+l -'"/tj H z h- 

1 



1 



Aa^c 



Ok 

4 



2 ' c 1 — a^c^ 
1 1 e'''+'+'''. 



Let 
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1 



1 



+ ■ 



Eliminating r[, we have 



2^2 



5yt 5yt-l 



+ V 



1 



1 



c 1 — a^c^ 
4a2c2 



+Oyt-l \ ; = 0, 



(3.32) 



c 1 — cP-c^ 



{Wk+l-Wk). 



This is a semi-discrete Camassa-Holm equation. Note that the lattice parameter depends on 



the time and space. 

Differentiating the first equation in (13.321) with respect to t, we obtain 

'Wk+l-Wk Wk-Wk-l 



□ 



-2 



+ 1 



4 



5^- 



dtWk+i+dtWk , 5 dtWk + dtWk-i 



Wk+i+Wk ^ 1 52 



+ 1 



5ti 



c 1 — a^c^ 

4a^c^ \ 



J 



1 + - 



Wk 



+ Wk-i , 1 



+ - 



c 1 — a^c^ 



V 



0, 



Taking a continuous limit, this leads to the CH equation. 
Then we have the following theorem. 

Corollary 3.4. The semi-discrete CH equation has a determinant form ofN-soliton solutions. 

Proof. From Theorem 3.2 and 3.3, the proof is obvious. □ 



4. Peakons, solitons and cuspons in the semi-discrete CH equation 

4.1. Peakon limit in the semi-discrete Camassa-Holm equation 

It is known that the CH equation has a peakon solution in the limit K ^ 0. Here we consider 
a peakon limit in the semi-discrete CH equation (|3.27l) . 
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In the semi-discrete CH equation (13.271 ). consider the limit 



5^ 



where 5(fc) is a dirac delta function. In this limit, the first equation of the semi-discrete CH 
equation (|3.27l) leads to 

d\w — w = d{X — cT) . 
This differential equation has a solution 

N 

w = J^A,-(r)exp|X-cr|, 

1=1 

which is a form of the peakon solution. 

Thus if we consider very small K with very small 5yt at some points k, the solutions of the 
semi-discrete CH equation tend to the peakon solutions of the CH equation. 



4.2. One soliton/cuspon solution 

From the determinant formula with a, i/a, 2 = ±1, the x-functions for one soliton/cuspon 
solution are 



(4.1) 



c + p J \c — p 

with = 2p{x — vt —xq), V = l/(c^ — p^) where c = 1/k^ > 0. This leads to a solution 

Ap^cv 

X = 2cJC + log(^|) , T = t, (4.3) 

where the positive case in Eq. (l4.2l) stands for one smooth soliton solution when p < c, while 
the negative case in Eq. (l4.2l) stands for one-cuspon solution when p > c. Otherwise, the 
solution is singular. Thus Eq. (|4.2|) for nonsingular cases can be expressed to 

Ap^cv 



^(^'^) (c^+p^) + \c^-p^\coshQ 
Similarly, for the semi-discrete case, we have 



c — p 



c + p 



l+ap\ 
\ —ap 



c + p 



c — p 



l+ap\ 
l—ap 



with = —2pv{t+xo), resulting in a solution of the form 

Ap^cv 



(c^-f- p^) + \c^ — p^\ 



\+ap 
l—ap 



-k 



(4.4) 



(4.5) 



(4.6) 



in conjunction with a transform between an equidistance mesh (a) and a non-equidistance 
mesh 

{\+ac)gk+ihk-{\-ac)gkhk+\ 



(1 +ac)gk+\hk+{\-ac)gkhk+\ ' 
Eq. (14.61) corresponds to 1-soliton solution when p < c, 1-cuspon solution when p > c. 



(4.7) 
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4.3. Two soliton/cuspon solutions 



From the determinant formula with ai^ijai^i = ±1, the T-functions for two soUton/cuspon 
solution are 



goc 


1 + 


ci -pi 
ci+pi 


e^' + 


ci-pi 

C2+P2 


e^^ + 


h 


1 + 


Cl +Pi 
ci -pi 


e'^ + 


C2+P2 
C2-P2 





(Ci -Pl){c2-P2) 



{ci+pi){c2 + P2) 

{ci+pi){c2 + P2) 



(Ci -Pl){c2-P2) 



P^-P2 \ ^01+02 
PI+P2J 

PI+P2J 



with 81 = 2pi{x-Vlt -Xio), 02 = 2p2{x-V2t -X20), Vl = l/{cj-pl), V2 = l/(c|-/7^). 

The parametric solution can be calculated through 



w(x,0= (log|)^, X = 2cx + log(£j, T = t, 



(4.8) 



whose form is complicated and is omitted here. Note that the above expression includes two- 
soliton solution (pi <c\,p2< 02), two-cuspon solution {p\ > c\,p2> 02), or soliton-cuspon 
solution {p\ <c\,p2> C2). 

Similarly, for the semi-discrete case, we have 



Cl -p\ 



c\+p\ 



1 +a/7i 
1 —ap\ 



.9' + 



C2-P2 



C2+P2 



+ 



{Ci-pi){c2-p2) 



{ci+pi){c2 + P2) 



PI -P2 
P1+P2 



n+ap2 y^Q, 

\l-ap2J 

l+api\ /'l+ap2\' 
l-apij \l-ap2J 



,61+62 



%ocl + 



ci+pi 



ci-p\ 



1 +api 
1 — api 



+ 



C2 + P2 



C2-P2 



l±ap2\ e, 



+ 



{C\+P\){C2+P2) 



{Ci-Pi){c2-P2) 



P\-P2 



P1+P2J \l-api 



l+api 



1 -ap2 
k 



^+aP2 ^ ^9i+62 
1 -ap2 



with 01 = 2pi {—vit—xio), 02 = 2p2{—V2t — X20) . The solution can be calculated through 



log — 

hkJt' 



Wk{t) -- 
with a transform 

{\+ac)gk+\hk+{\-ac)gkhk+\ ' 
The form is complicated and is omitted here. 



{l+ac)gk+ihk-{l- ac)gkhk+i 



(4.9) 



(4.10) 



5. Numerical computations 

In this section, several examples will be illustrated to show the integrable semi-discretization 
of the CH equation is a powerful scheme for the numerical solutions of the CH equation. They 
include (1) Propagation of one-cuspon solution; (2) Interaction of two-cuspon solutions; (3) 
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Head-on collision of soliton-cuspon. In actual computations, the initial mesh spacing is 
assigned by 

(l+ac)4)^+i-(l-ac)4>fe 

^'^ = 2-— r- — ; r— , (5.1) 

[l+ac)^k+i + {}-ac)^k 
where (|)^ = gk/h]^ is obtainable from the corresponding determinant solutions. At each time 
step, from the second equation in (13.271) . the evolution of can be exactly calculated by 



dl+'^i:^, , (5.2) 

with = (2 + 6^)/(2 — 5^). Then, the solution w^^^ is easily computed by solving a 
tridiagonal linear system based on the first equation of the scheme. Therefore, the computation 
cost is less than other traditional numerical methods. Furthermore, we comment that the 
mesh spacing 5^^ is driven automatically by the solution during the numerical computation. 
Therefore, we would like to call it the self-adaptive method. 

Example 1: 1-cuspon propagation. The parameters taken for the 1-cuspon solution 
are p = 10.98, c = 10.0. The number of grid is taken as 100 in an interval of width of 4 
in the .)c-domain, which implies a mesh size ofh = 0.04. However, through the hodograph 
transformation, this corresponds to an interval of width 74.34 in theX-domain, which implies 
an average mesh size of 0.7434. In subsequent examples except for Example 4, the grid 
number is fixed to be 100. The time step size is taken as At = 0.0004. Fig. 1(a) shows the 
initial condition. Figs. 1 (b)-(c) display the numerical solutions (solid line) and exact solutions 
(dotted line) at t = 2,4, respectively. The norm are 0.0365 at f = 2, and 0.0985 at ? = 4. 
It is noted that the numerical error is mainly due to the numerical dispersion. In other words, 
even after a fairly long time, the numerical solution of a cuspon preserves its shape very well 
except for a phase shift. The reason for the numerical dispersion is thought due to the explicit 
method to solve the equation for the time evolution of mesh size. The detailed analysis is left 
for our another paper focusing on numerical solutions of the CH equation. 

Example 2: 2-cuspon interaction. The parameters taken for the two-cuspon solution 
are pi = 11.0, p2 = 10.5, c = 10.0. Fig. 2(a) shows the initial condition, and Figs. 2(b)- 
2(e) display the process of collision at several different times. As far as we know, what is 
shown here is the first numerical demonstration for the cuspon-cuspon interaction due to the 
singularities of cuspon solutions. As shown in Fig. 2(e), the 2-cuspon solution regain their 
shapes after the collision, only resulting in a phase shift. As mentioned in [18], the two cuspon 
points are always present during the collision. 

Example 3: Soliton-cuspon interaction. Here we show two examples for the soliton- 
cuspon interaction with c = 10.0. In Fig. 3, we plot the interaction process between a soliton 
of pi = 9.12 and a cuspon of p2 = 10.98 at several different times where the soliton and 
the cuspon have almost the same amplitude. It can be seen that another singularity point 
with infinite derivative (Wx) occurs when the collision starts (t = 12.0). As collision goes 
on (t = 14.4, 14.6, 14.8), the soliton seems 'eats up' the cuspon, and the profile looks like a 
complete elevation. However, the cuspon point is present at all times, especially, att = 14.6, 
the profile becomes one symmetrical hump with a cuspon point in the middle of the hump. 
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Figure 1. Numerical solution of one single cuspon solution: (a) t — 0.0; (b) t — 2.0; (c) t — 4.0. 



In Figure 4, we present another example of a collision between a soliton {p\ =9.12) and 
a cuspon {p2 = 10.5) where the cuspon has a larger amplitude (2.0) than the soliton (1.0). 
Again, when the collision starts, another singularity point appears. As collision goes on, the 
soliton is gradually absorbed by the cuspon. Att = 10.3, the whole profile looks like a single 
cuspon when the soliton is completely absorbed. Later on, the soliton emerges from the right 
until t = 16, the soliton and cuspon recover their original shapes except for a phase shift when 
the collision is complete. 

Example 4: Initial condition of non exact soliton solutions. We show that the 
integrable scheme can be also applied for the initial value problem starting with a non- 
exact solution. We choose an initial condition in the following procedure. The mesh size 
is determined by 

5^ = 2ac(l - 0.8sech(2yta - Wjc/3)) , (5.3) 

where (= 8) is the width of computation, N = 201 is the number of grid in x-domain, 
^=1,---,A^— 1, and a = W^/ {N—l)= 0.04. Then the initial profile can be calculated through 
the second equation of the semi-discretization (13.271) . The initial profile is plotted in Fig. 5(a). 
Figs. 5(b), (c) and (d) show the evolution at t = 10, 20, 30, respectively. It can be seen that 
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Figure 2. Numerical solution for the collision of two-cuspon solution with pi = 11 
P2 = 10.5, c = 10.0: (a) t = 0.0; (b) t = 13.0; (c) t = 14.8; (d) t = 16.6; (e) t = 25.0. 
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Figure 3. Numerical solution for soliton-cuspon collision with p\ 
c = 10.0: (a) t = 0.0; (b) t = 12.0; (c) t = 14.4; (d) t = 14.6; (e) t 
t = 25.0; 



= 9.12, p2 = 10.98 and 
= 14.8; (f) t = 17.0; (g) 
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Figure 4. Numerical solution for soliton-cuspon collision with pi = 9.12, p2 = 10.5 and 
c = 10.0: (a) t = 0.0; (b) t = 9.0; (c) t = 10.0; (d) t = 10.3; (e) t = 10.6; (f) t = 11.5; (g) 
t = 16.0. 
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Figure 5. Numerical solution starting from an initial condition: (a) t = 0.0; (b) t = 10.0; (c) 
t = 20.0; (d) t = 30.0. 

a soliton with large amplitude is developed first, and moving fast to the right. By t = 30, a 
second soliton with small amplitude is developed and a third soliton is bom from the second 
soliton. 

6. Concluding remarks 

An integrable semi-discretization of the CH equation has been presented in this paper. 
Determinant formulas of the A^-soliton solutions of both the continuous and semi-discrete CH 
equations have been derived. Multi-soliton, multi-cuspon and multi-soliton-cuspon solutions 
can be generated from above determinant formulas. As further topics, we attempt to construct 
integrable semi-discretizations of the Degasperis-Procesi equation and other soliton equations 
possessing non-smooth solutions such as peakon, cuspon, or loop-soliton solutions. We will 
address these issues in forthcoming papers. 

Applying integrable discretizations of soliton equations to numerical computations 
remains a promising but not thoroughly explored subject. In the present paper, even for 
relatively large mesh sizes, very accurate numerical solutions of cuspon-cuspon and soliton- 
cuspon interactions for the CH equations are achieved through our proposed integrable semi- 
discrete scheme. In addition, a numerical computation starting with an non-exact initial 
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condition is performed with a satisfactory result. It is worth noting that the integrable semi- 
discrete scheme of the CH equation is also a self-adaptive method, which is of great interest 
in the area of numerical partial differential equations. We intend to extend this new idea of 
self-adaptive method to other PDEs in the near future. 



7. Appendix 



In this appendix, we prove the differential formulas for x„, (|2.9I )- (|2.16I) . Let us introduce a 
simplified notation. 



(7.1) 



In this notation, we have x„ = |\|/("),\|/("+^), ■ ■ ■ ,\[;("+^ ^)|, thus differentiating x„ by x and 
using (|2.3I) we obtain 

axx„ = , v^(''+i) , ■ ■ ■ , a.v^"+^'^ , ■ ■ ■ , ^^"^"-'^ I 

N-l 

= £ |\|/W,V"+i),...,V"+^'+i) + c\|/("+^'\ 
= |v|/W, • • • , ■ ■ ■M"^''-'^ I 



(n+W-l)! 



A'-l 



(n) ,,f("+l) ... ,,f(«+A'-l) 



which gives (12.91) . Similarly differentiating x„ by j and using (12.41) we get 



A'-l 



A'-l 



+ 2c|V"W"+^\ ■ ■ ■ , V"+^"2)^^(n+A')| ^^^2|^(n)^^(n+l)^ . . . ^ ^(n+A^-l) | _ 

We obtain (12.101) from the above equation and (12.91) . Eqs. (12.111) and (12.121) can be proved 
by using (12.51) and (12.61) in a similar way. Finally eqs. (I2.13l) - (|2.16l) can be verified by 
differentiating (12.91) and (|2.10|) by t and through similar calculations. 
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